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The distribution of thermally excited vortex loops near a superfluid phase transition is calculated 
from a renormalized theory. The number density of loops with a given perimeter is found to change 
from exponential decay with increasing perimeter to algebraic decay as T c is approached, in agree- 
ment with recent simulations of both cosmic strings and high-T c superconductors. Predictions of 
the value of the exponent of the algebraic decay at T c and of critical behavior in the vortex density 
are confirmed by the simulations, giving strong support to the vortex-folding model proposed by 
Shenoy. 
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The role of thermally excited vortex loops in three- 
dimensional phase transitions where a U(l) symmetry is 
broken recently has become a prime topic in cosmology 
0U and high-T c superconductivity |||-||]. These tran- 
sitions arc in the same universality class as the super- 
fluid A-transition in 4 He. In the helium case our orig- 
inal renormalization theory based on vortex loops ^ 1 } 
has been extended to calculations of the specific heat S ] 
and to the dynamics of the transition fl. In this theory 
the Landau- Ginzburg- Wilson Hamiltonian is rewritten to 
cast it in terms of its elementary excitations, spin waves 
and vortex loops, providing an alternative method for 
carrying out the renormalization process, compared to 
the more traditional perturbation theories which expand 
the Hamiltonian expressed in terms of the order param- 
eter. 

Here we further employ the loop theory to gain insight 
into recent simulations of the high-T c transition and 
of cosmic-string phase transitions in the early universe 
0,§. The probability of occurrence of a vortex loop with 
perimeter P is calculated, and in agreement with the 
simulations we find a crossover from quasi-exponential 
decay of the probability with perimeter at low temper- 
atures, to purely algebraic decay precisely at T c . This 
provides strong support for the phcnomcnological " Flory 
scaling" treatment of the random-walking loops devel- 
oped by Shenoy and co-workers |Tc|| . 

In the vortex-loop theory the superfluid density p s is 
reduced by thermally excited loops whose average diam- 
eter a increases as the temperature is increased, and the 
density is finally driven to zero at T c by loops of infinite 
size 0j. Defining a dimensionless superfluid density by 
K r = h 2 p s a /m 2 kBT, where m is the mass of the 4 He 
atom and a a the smallest ring diameter, the equation for 
the renormalized density is given by (§] 



Here A Q = 47r 3 /3, K d is the "bare" superfluid density 
resulting from the spin waves (and is the initial value of 
K r ), and U(a) is the renormalized energy of a ring, given 
by 
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where C is a nonuniversal constant characterizing the 
core energy. For helium C and a Q are determined from 
two experimental inputs, T c = 2.172 K and the amplitude 
of the superfluid density yielding C = 1.03 and a D — 
2.5 A. The effective core size a c in Eq. (2) was suggested 
by Shenoy and co-workers @,|l(J to be a result of the ran- 
dom walk of the loop giving rise to radial fluctuations of 
order a c about the average diameter. This folding of the 
loop occurs because antiparallel vortex segments lower 
the energy. A simple polymer-type calculation fi"c[ | us- 
ing energy-entropy arguments yields a c /a — (K r a/a ) e , 
where 8 = d/(d + 2) = 0.6 in d = 3 dimensions has the 
same form as the well-known Flory exponent of the self- 
avoiding walk. 

Eqs. (1) and (2) then constitute a coupled set of in- 
tegral equations for the renormalized superfluid density, 
and can be solved recursively starting from the bare scale 
a and iterating to distances greater than the correlation 
length £ = a /K r . In practice these are converted to a set 
of coupled differential equations similar to the Kosterlitz 
recursion relations pd[ | for the two-dimensional case, and 
are solved using a Runge-Kutta technique ||. As T is 
increased (K Q decreased) the solution for p s falls to zero 
as (T c -T) v , with v = 0.6717 for 6 = 0.6. This can be 
better matched to the most precise experimental value 
[pf v = 0.6705 by adjusting to 6 = 0.594, which is rea- 
sonable since it is known that the Flory-type arguments 
are not exact in three dimensions [fl3| . 

The arguments of Ref. 10 also yield a result for the 
average perimeter of a loop of diameter a, 
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where B is a constant and the exponent S — 2/(D+2) 
= 0.4. This form for the perimeter was at least partially 
verified in the computer simulations of Ref. 10, and in the 
XY model simulations of Epiney jl4| . For P/a a greater 
than about 20 the Epiney data for the average loop size 
versus average perimeter can be fit by Eq. (3) with B = 
1.8, although the resolution is poor because of scatter in 
the data resulting from the relatively small lattice (16 3 ) 
that was simulated. 

The distribution of the density of loops with a given 
perimeter P can be obtained from the theory outlined 
above, which is of interest because these distributions 
have now been measured in the cosmic-string and 
high-T c simulations. The probability per unit vol- 
ume for finding a loop of mean diameter between a and 
a+da is 
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Equating this to the probability D(P) dP/a of find- 
ing the corresponding loop of perimeter between P and 
P+dP gives the probability distribution 
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For temperatures well below T c , U(a) ~ a In a, and hence 
D(P) decreases exponentially with P . Near T c , how- 
ever, the behavior is quite different. By differentiating 
Eq. (1) with respect to a and substituting for the expo- 
nential term in Eq. (5) gives 



D(P) 



2alB 



d (1/K r ) 
da 



(6) 



Precisely at T c , K r from Eq. (1) satisfies the condition 
K r a/a = D Q = 0.39, at least for values of a greater than 
about 5a , and where D Q is a universal constant |p|,|l5|. 
This condition is the three-dimensional equivalent of the 
univeral "jump" of the superfluid density in two dimen- 
sions fl"l"|| . Inserting this result into Eq. (6) and employing 
Eq. (3) yields the prediction that at T c the loop distribu- 
tion will cross over from exponential to algebraic decay 
with P , 
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where the exponent 7 = 3£ + 1. For the "Flory" value 
5 = 0.4 this would predict 7 = 2.2. This form for D(P) 
signals the onset of loops of infinite size, since they no 
longer have a vanishing probability. It should be noted 
that the algebraic decay is a consequence of the strong 



renormalization at T c , where the screening of large loops 
by smaller ones causes the variation of U(a) to change 
H from ~ a In a to In a at T c , causing the change from 
Eq. (5) to Eq. (7). 

The crossover from exponential to algebraic decay is 
a central feature observed in the recent cosmic-string 
and high-T c simulations using lattices larger than 
96 3 , and which was also seen with more limited resolution 
in the earlier XY model simulations of Epiney jl4j . The 
results of Antunes and Bettencourt [|| yield 7 = 2.23 
± 0.04 at T c , which from the analysis above gives S = 
0.41 ± 0.01, in very good agreement with Shenoy's Flory- 
scaling prediction. Fits to the high-T c simulations of 
Nguyen and Sudbo (in the zero-field, isotropic limit 
of their Villain model) give 7 = 2.4 [M, leading to a 
higher value 5 — 0.48. However, the results of Ref. 2 
show that 7 increases rapidly above T c to the value of 
2.5 found by Vachaspati and Vilenkin fl7j| , and hence an 
accurate determination requires bracketing temperatures 
very close to T c . It is interesting that the result 7 — 
3 5 + 1 apparently remains valid even above T c , since 
the Vachaspati- Vilenkin value of 7 = 2.5 is based on the 
"Brownian" value 5 = 0.5. 

The magnitude of the loop distribution D(P) at T c as 
calculated from either Eq. (5) or (7) appears to be about 
a factor of 3 smaller than found in the simulations. Com- 
paring the continuum calculation to the lattice results is 
made difficult by uncertainties in matching at the scale 
of the lattice spacing ai; in computing the magnitude of 
Eq. (7) for the comparison the geometric value a = \[2 ai 
was employed, but it is not entirely clear that this is the 
correct choice. 

The total length per unit volume p v of the vortex loops 
can be found by multiplying D(P) by P and integrating. 
At T c this can be found explicitly using Eq. (7), 



Pv(T c 



ttSB 



2a 2 A D {3S-l) 



(8) 



The quantity p v (T c )af was postulated in Ref. 1 to be 
universal, with a value of 0.6 in lattice units (0.2 per 
placquette, with three placquettes per unit volume). This 
means that the coefficient B in Eq. (8) characterizing 
the relationship between the perimeter and the average 
loop diameter must be universal, since all of the other 
parameters are. As with D(P) above, the magnitude of 
Eq. (8) must be multiplied by a factor of about 3 to 
match with the lattice results. For liquid helium, the 



vortex density |18| at T c is then about 1 / a„ ~ 1 x 10 15 



/cm 2 , which is many orders of magnitude higher than 
previous estimates Jl9| which did not use a renormalized 
theory. 
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FIG. 1. Normalized vortex density as a function of T c / T, 
for several different lattice sizes. 

A further prediction of the loop theory is the exis- 
tence of critical behavior in the vortex density at T c [[l8| , 
which has now been seen in the cosmic-string simula- 
tions This is shown in Fig. 1, which plots the nor- 
malized density versus T c / T, calculated by integrating 
Eq. (5) very near the transition. The behavior above T c 
(dashed lines) is only conjectured, as the theory is not 
valid there. To suppress the exponential variation aris- 
ing from thermal excitation of the smallest loops, the core 
energy constant was taken to be the relatively large value 
C = 4/3, corresponding to the Villain model. The den- 
sity just below T c is found to decrease from its value at 
T c as -{T c -T) s . This exponent was not measured in the 
cosmic-string simulations, probably because of the rela- 
tively low core energy in that model, which would make 
it difficult to separate the algebraic behavior from the 
exponential. However, the exponent of the density in the 
region just above T c was measured in Ref. 1 to be 0.39 
± 0.01, which is quite consistent with the value of 8 = 
0.41 determined above from the same type of simulation. 
It is well known from scaling and renormalization-group 
studies that many critical exponents are identical above 
and below the transition, and it is likely that this expo- 
nent constitutes a further measurement of 8 . Since 8 is 
less than one, T c marks an inflection point in the density, 
where the curvature changes sign. 

Also shown in Fig. 1 is the effect of a finite-size lattice 
on the density near T c . For this the recursion relations 
are stopped at a finite scale that is a fraction j3 — 0.75 
smaller than the lattice size, where /3 was found in Ref. 8 
by comparing to the helium simulations of Ref. 15. The 
effect of finite size is to smear out the critical behavior, 
and this explains why it was not apparent in the early 
simulations pfj| , which only used a maximum 10 3 lattice 
size. 

The coherence length in the loop theory can be iden- 
tified with the diameter of the largest loop that is ther- 
mally excited flq] . Since this is divergent at T c , in a 



system of finite (but macroscopic) size the transition can 
be identified with the point where a loop just touches 
two opposing boundaries of the system. This is the per- 
colation threshold for establishment of the infinite vortex 
cluster, as first proposed by Onsager 50 years ago, and 
verified in simulations [ pT| . At higher temperature even 
larger loops can be excited; these can take the form of a 
single vortex line passing from one side of the system to 
the other, with the topological return path of the loop be- 
ing around the outside of the system. These are known as 
"string" excitations in cosmology p]^,p7|. The left-over 
strings from the rapid cooling through the phase transi- 
tions in the early universe are thought to be the source of 
matter |l|,|2^] . A laboratory example of this may be the 
observation of a few remnant vortex lines (23) in a finite- 
sized sample of liquid 4 He which has not been rotated or 
stirred, but simply cooled through the A-transition. 

The key role of thermally excited vorticity in the 
present model of the phase transition allows a rather dif- 
ferent viewpoint of the Kibble-Zurek mechanism |[9 22 
for defect formation in rapidly quenched transitions as 
discussed above. This is of interest because of the possi- 
bility of carrying out such experiments in liquid helium 
]l^,||j24U2|]. In this view the vortices "created" in a 
quench of superfluid 4 He through T c are simply a per- 
turbation on the equilibrium vortex density. This pertur- 
bation is a consequence of the dynamics of large loops, 
which become slow in their response to external fields, 
and which are actually the source of the critical slowing 
down of the equilibrium transition || . In a rapid quench 
the large loops in the vicinity of T c cannot keep up and 
fall out of equilibrium, forming an excess density that 
survives to lower temperatures, finally decaying to the 
equilibrium line density only after the quench is stopped. 
It may be possible to model this process analytically us- 
ing the Fokker-Planck equation for the loop distribution 
function formulated in the dynamic theory of Ref. 9, with 
the temperature being a ramp function in time. It should 
be noted that the exponent 8 from above plays an impor- 
tant role in the vortex dynamics, since the frictional drag 
force on a loop is proportional to its total perimeter. It 
was found in Ref. 9 that the dynamic exponent charac- 
terizing critical slowing down is given by x = z where z 
= {1-8)/ 8. For the "Flory" value 8 = 2/(d+2) this gives 
the scaling result z = d/2 = 3/2; for <5 = 0.41 as found 
above, z = 1.44, a few percent smaller. The possibility 
of deviations from the scaling result has been suggested 
previously in perturbative dynamic theories p6| . 

Vortex creation has been observed in superfluid 3 He 
p5| where the quench is induced by absorbed neutrons 
depositing their energy in a small region of the liquid, 
which heats it above T c , and which is then rapidly cooled 
back down by the surrounding cold liquid p5[ . Although 
3 He is a p-wave BCS superfluid, vortex loops will still be 
associated with the phase transition as above, but due to 



the Ginzburg criterion |27| they will only be appreciable 
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in a very narrow temperature range near T c , since the 
zero-temperature coherence length of 3 He is of order 500 
A, compared with a = 2.5 A in 4 He. The theory p3] of 
the quenched 3 He involves the exponents <5 and 7 at T c ; 
the use of the Brownian values [[l7| 7 = 2.5 and S =0.5 
needs to be reexamined in light of the present results. 

The vortex loop model also allows insights into the 
high-T c superconducting transition in zero field. If T c 
is the point where the loops of infinite size act to bring 
all supercurrents to a halt, then it is not to be identi- 
fied as the point where thermal de-excitation of Cooper 
pairs is complete. The continued existence of pairs above 
T c has been suggested in experiments, commonly known 
as the pseudogap phenomenon The vortex loops 

constitute a concrete physical picture of the " phase fluc- 
tuations" postulated to give rise to this effect in Ref. 28. 
The Cooper pairs above T c will not be the same as those 
below, since they will no longer be part of a macroscopic 
BCS-type condensate, which is destroyed by the vortices. 
Presumably the pairs would be more localized excita- 
tions, on the scale of the 10-15 A zero-temperature co- 
herence length. 

The vortices also offer a simple explanation |2^] of the 
magnetic flux noise in YBCO samples that is observed 
to increase rapidly by many orders of magnitude over a 
temperature range of order 5 K near T c |30| . When the 
loops being thermally excited terminate on the sample 
surface, they will induce dipolar current patterns on the 
surface, and this will generate magnetic flux that can 
be sensed by a detection loop above the surface. As T c 
is approached from below both the number and size of 
the loops will increase, increasing the flux noise through 
the detector. This effect can also be observed in a low-T c 
superconductor, since the same vortex-loop transition oc- 
curs also in that case, but with the considerable difference 
that the large zero-temperature coherence length (several 
thousand A) causes the temperature range where the vor- 
tices are appreciable to be very much closer to T c . The 
experiments |3l| showed that indeed a very sharp flux- 
noise peak could be observed in a lead sample only within 
about 2 mK of T c , and that this was only an upper limit 
to the width due to the resolution of the thermometry 
and the additional broadening that would be caused by 
a distribution of T c 's across the sample. 

In summary, a vortex-loop theory is able to provide 
physical insight into recent models of cosmic strings and 
high-T c superconductors. The theory relates the critical 
exponents measured in the simulations to the random- 
walk nature of the loops, and the good agreement be- 
tween the predicted and measured exponents provides 
strong support for the Flory-scaling ansatz of Shcnoy . 
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